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Abstract 

We have calculated the optical absorption for InGaNAs and GaNSb using the band anticrossing (BAG) model and a 
self-consistent Green's function (SGGF) method. In the BAG model, we include the interaction of isolated and pair N 
levels with the host matrix conduction and valence bands. In the SGGF approach, we include a full distribution of N 
states, with non-parabolic conduction and light-hole bands, and parabolic heavy-hole and spin-split-off bands. The 
comparison with experiments shows that the first model accounts for many features of the absorption spectrum in 
InGaNAs; including the full distribution of N states improves this agreement. Our calculated absorption spectra for 
GaNSb alloys predict the band edges correctly but show more features than are seen experimentally. This suggests 
the presence of more disorder in GaNSb alloys in comparison with InGaNAs. 

Keywords: Dilute nitride semiconductors; Self-consistent Green's function; Optical absorption; 
Band anticrossing model 



Background 

The substitution of a small fraction x of nitrogen atoms, 
for group V elements in conventional III-V semicon- 
ductors such as GaAs and GaSb, strongly perturbs the 
conduction band (CB) of the host semiconductor. The 
band structure of dilute nitride alloys has been widely 
investigated [1]. We have recently developed [2] a SCGF 
approach to calculate the density of states (DOS) near the 
conduction band edge (CBE) in these alloys. 

One way to test the accuracy of this model is to look 
at optical absorption spectra for dilute nitride samples, 
where we expect to see features related to the N states 
present in the samples. The absorption spectrum arises 
from transitions between valence and conduction band 
states. It provides knowledge of the energy gap in semi- 
conductors, and also gives significant information about 
the band structure of materials. Experimental measure- 
ments of absorption spectra can be used to benchmark 
band structure calculations. In this paper, we investigate 
two different materials: In^Gai.yN^^Asi.^c, for which the 
band structure has been widely studied and many of the 
features are well established and GaN;cSbi_;c> for which 
much less information has been reported in the literature. 
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We consider two different models for the band struc- 
ture of dilute nitride alloys, firstly a five-level band anti- 
crossing (BAG) model, including the host semiconductor 
CB and valence bands, isolated N and pair N-N states 
and, secondly, the linear combination of isolated nitrogen 
states (LCINS) model [3,4], which allows for interaction 
between N states on nearby sites as well as inhomoge- 
neous broadening and produces a distribution of N state 
energies. In the LCINS model, the band structure of the 
alloys is calculated using a SCGF approach [2]. 

For In^Gai-^/N^cAsi-^c alloys, we find that the BAG model 
reproduces the main features in the absorption spectrum, 
in agreement with previous work [5,6]. However this 
model shows some additional features which are related to 
the N and N-N state energies, reflecting that in the BAG 
model, we have ignored the actual distribution of localised 
states. Including the LCINS distribution of N states in 
InyG^Li-yNxAsi-x using the SCGF approach [2] removes 
the additional features found in the BAG calculations and 
gives absorption spectra that are in very good agreement 
with experimental data. 

We then apply our methods to GaN;cSbi_;c> where much 
less information is known theoretically and experimen- 
tally. The overall width of the optical spectrum can be well 
fitted by our models for the absorption spectrum. Both the 
BAG and LCINS models account for the absorption edge 
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of GaNSb alloys, supporting the presence of a band anti- 
crossing interaction in these alloys. However, the five-level 
BAG model gives more features than are seen experimen- 
tally in the absorption spectrum. Including a distribution 
of localised state energies, obtained by modifying those 
calculated for GaNAs, makes the calculated absorption 
spectra smoother and gives better agreement with exper- 
imental data but still shows some discrepancies around 
the localised state peak energies. These results suggest 
the presence of more disorder in GaNSb samples than in 
InGaNAs. This disorder may be due to sample inhomo- 
geneities or due to an intrinsically broader distribution of 
N states in GaNSb than in InGaNAs. 

The remainder of this paper is organised as follows. In 
the 'Methods' Section we first provide an overview of opti- 
cal absorption calculation, followed by a description of the 
band structure models used for dilute nitride alloys. 
The theoretical results for InGaNAs and GaNSb are pre- 
sented and compared with experiment in the 'Results and 
discussion' Section. Finally, we summarise our conclu- 
sions in the last section. 

Methods 

The absorption spectrum a(E) describes the rate of 
absorption of photons with energy E = hco per unit 
distance and can be described using a single-electron 
approximation. In this approach, the absorption spec- 
trum for allowed transitions between valence band v and 
conduction band c states is given by [7-10] 



acvihco) = 



Y,MlS{Eck-E,k-hxa), 



(1) 



where fm is the photon energy, e and mo are the elec- 
tron charge and mass, c is the speed of light and rir 
is the refractive index. / is the Fermi-Dirac distribution 
function. Here, we assume a filled valence and empty con- 
duction band, so that/v — / = 1. The matrix element in 
Equation 1 can be written for transitions between valence 
p- and conduction 5-like zone centre states as [9] 



Mi 



moEp 



(2) 



where Ep = (2mo/ff)p^ is the interaction energy, and the 
momentum interband matrix element,/?, can be estimated 
from experiment as 



^ mo) 2m* (Eg + I Aso)' 



(3) 



where m* is the conduction band effective mass. Eg is the 
band gap (between conduction and valence band) and Ago 



is the spin-orbit-splitting energy. Therefore, Equation 1 
can be written as 



Tte^h Ml . 



(4) 



where E = E^k — Eyj^ is the transition energy between con- 
duction {Eck) and valence {Ey^) states with wavevector ky 
Jcv{E) is the joint density of states, and we ignore for now 
the energy dependence of M^. 

The hole-electron interaction can be included assuming 
Elliot's theory, which applies to parabolic and nondegener- 
ate bands. According to Elliot's model [11], the absorption 
spectrum is modified because of the hole-electron inter- 
action through a multiplicative function F^x given by [8] 



271 g 



1 — exp(— 2;r^) 



(5) 



where 



\h(j) — 



1/2 



(6) 



and Ry is the exciton Rydberg energy, given by [12,13] 



-^(^) = 13.61 (^)eV, (7) 



where k is the static dielectric constant, h is the Planck 
constant, /x = l/{m~^ + m~^) is the reduced mass, and 
niy = {n^il^ + m^^^)^/^, where mpi and m/ are the heavy- 
hole and light-hole effective masses. 

Including the hole-electron interaction, the total 
absorption spectrum atot(E) can then be obtained as 



(8) 



where of//^, of/^/^ and aso are the absorption spectra for tran- 
sitions from the light-hole (LH), heavy-hole (HH) and 
spin-orbit split-off (SO) bands to the conduction band, 
respectively. 

The effect of the incorporation of N in (In)GaNAs alloys 
can be described in different ways. We investigate here 
how the model chosen influences the calculated alloy 
absorption spectrum. We first present a simple model 
including isolated and pair N states using the BAG model. 
This model includes the nonparabolicity of the conduc- 
tion and light-hole and split-off bands and the interaction 
between the split-off and conduction bands. In the sec- 
ond model, we then include the full LCINS distribution of 
localised states using the SCGF model. 

Optical absorption of dilute nitride alloys in five-level BAG 
model 

Here we first consider a simpler model, including iso- 
lated and pair N states and their interaction with the host 
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semiconductor conduction, valence and spin-orbit split- 
off bands. The conventional BAG model treats the host 
III-V conduction band dispersion as a parabolic band [14]. 
Test calculations that we have undertaken show that the 
inclusion of band non-parabolicity can strongly modify 
the calculated absorption spectra due to the change in 
the joint density of states caused by the non-parabolicity. 
In order to treat the host matrix conduction band non- 
parabolicity and the effect of N on the alloy conduction 
band dispersion, we construct a 5 x 5 Hamiltonian. This 
includes the Kane nonparabolicity of the host matrix con- 
duction band, due to interactions with the light-hole and 
split-off bands, and treats the effects of N using a three- 
level BAG model [15], including isolated and pair N states. 
This Hamiltonian is given as 



Table 1 BAC model parameters for InyGai-yNxAsi-x 







Vnc Vnnc 










Vnc 


En 0 


0 


0 




H = 


Vnnc 


0 Enn 


0 


0 


(9) 






0 0 




0 








0 0 


0 EyQ 


- ^so 





where A^^ indicates the spin-orbit-splitting energy, and 
EyQ is the energy of the valence band maximum. The 
energy of the isolated N levels {E^i), N-N pair states {E^n) 
and the conduction band edge (£co) are assumed to vary 
with composition, Xy and temperature, T, as [15] 



Enn = Emno + Ynn^ + ^nnT, 



and 



(10) 

(11) 

(12) 



The interaction parameters are assumed to vary as Vjsic = 
Pns/^ and Vnnc = Pnn^/^nn> where the concentrations 
of single N and N-N pair states xm and xmn^ respec- 
tively, are determined from the total N concentration x as 

^A^A^ = and 

xn — ^ 2xmN' The chosen values of the 
above parameters for (In)GaNAs and GaNSb are given in 
Tables 1 and 2, respectively [16-18]. 

Galculating the eigenvalues of the matrix in Equation 
9 gives the dispersion for five bands, namely, the light- 
hole and split-off valence bands and three conduction 
bands denoted by £/, Em and Eu for lower, middle and 
upper bands, respectively. However, the five-band model 
of Equation 9 overestimates the LH nonparabolicity and 
omits the heavy-hole band dispersion. Figure 1 displays 
the band dispersion for Ino.04Gao.96No.01Aso.99 where we 



Parameter 


Symbol 


Values 


N energy 


Eno 


1 .706(1 - y) + 1 .44y - 0.38y(1 - y) (eV) 


N-N energy 


Enno 


1 .486(1 - y) + 1 .44y - 0.38y(1 - y) (eV) 


dEN/dT 


On 


-2.5 X 10-4 (eV/K) 


dENN/dT 


onn 


-2.5 X 10-4 (eV/K) 


dEN/dx 


YN 


-0.22 (eV) 


dENN/dx 


Ynn 


-0.22 (eV) 


dEc/dx 


Yx 


-2.1 eV 


N interaction 




1.97(1 -y)-F2y-3.5y(1 -y) (eV) 


N-N interaction 




2.69(1 -y)-F2y-3.5y(1 -y) (eV) 


Energy gap 


EcO 


%GaAs-1.33y-F0.27y2 (eV) 



have included a parabolic heavy-hole band, £hh = ^vO — 
h^k^/ (2m fi) and a nonparabolic light-hole band calculated 
using the six-band Luttinger-Kohn (LK) valence band 
Hamiltonian [21]. In this model, the band dispersion of 
the light-hole band is given by 



£/+£so , 
ElH = \- 



+ 3 



(13) 



where Ei = Eyo - h^k'^/(2mi) and Eso = Eyo - Ago - 
h^k^ /(2mso) are the parabolic light-hole and spin-orbit 
bands, respectively, and 



72 



mo 



(---) 

\mi mhj 



(14) 



The fractional F character, defined as the contribution 
of the host matrix GB states to a given level, can be cal- 

Table 2 GaNxSbi-x parameters at room temperature 

Parameter Symbol Values 



Lattice constant [19] 


oo 


6.09593 (A) 


Electron effective mass [1 9] 






Conduction 


m* 


0.039(mo) 


Liglit liole 


mi 


0.0439(mo) 


Heavy liole 


rrih 


0.25 (mo) 


Split-off 


mso 


0.12(mo) 


SO splitting energy 


Aso 


0.76 (eV) 


Energy gap [20] 


Eg 


0.725 (eV) 


Refractive index [19] 


rir 


3.8 


N energy [17] 


Eno 


0.82 - 2.3 X (eV) 


N-N energy [1 7] 


Enno 


0.48 - 2.3 X (eV) 


N interaction [1 7] 




2.4 (eV) 


N-N interaction [17] 




3.39 (eV) 
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Figure 1 The band dispersion for lno.04Gao.96No.01 Aso.99 in five-level BAC model. Arrows show the possible optical transition from spin-orbit 
split-off band (solid arrows), light hole (dashed arrows) and heavy hole (dash-dotted arrows) to conduction sub-bands. 



culated by finding the eigenvectors of Equation 9 and is 
given by 



fr(E) = 



1 + 



Vi 



Nc 



: + - 



V, 



NNc 



. (15) 



Here, for each conduction sub-band, we use the appro- 
priate energy £, as shown in Figure 1 by £/, £^ and 

The joint density of states for transitions from the 
valence band vi to the conduction band cf, times the host 
matrix F character of the conduction state, determines the 
absorption strength. It is given by 



Jcf,vi — 



7r2 (d£/,,/d^)' 



(16) 



where Ej^i = hco = Eyi — E^f is the energy separation 
between the CB cf and the valence band vi states. The 
absorption spectrum for transitions between bands vi and 
cf is given, using Equation 4, by 



Tze^h Ml 
acf,vi(hco) = :j^Jcf,yi(hco), 



(17) 



Figure 1 represents the band dispersion and all possible 
optical transitions for Ino.04Gao.96No.01Aso.99. The total 
absorption spectrum is given by the summation of nine 
individual absorption components, shown in Figure 1, for 
transitions from VBs vi (including LH, HH and SO bands) 
to CBs cf calculated using the 5-level BAC model {i.e. Ei, 
Em and Eu) as 

+ 0(u,LH + 0^/,SO + 0(m,SO + 0(u,So) > (18) 

where oiso-h c^so-m and aso-w are the absorption spectra 
from split-off band to lower, middle and upper sub-bands, 
respectively, and with a similar notation used for transi- 
tions from the HH and LH bands. 

Optical absorption of dilute nitride alloys in the LCINS 
model 

In order to include the full distribution of N states, we 
need to use the Greens function method in the frame- 
work of the LCINS model. In our recent work [2], we 
showed that the conduction band Greens function could 
be written as 



1 



E 



\Vj\^ 



E-Ej- AEj{E) 



(19) 
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where the (complex) energy shift of each locaUsed state / 
is given by 



(20) 



We can solve Equations 19 and 20 self-consistently by an 
iterative method [2] to calculate AEj(E). The density of 
states per unit volume for the CB, projected onto the host 
matrix conduction states, is given by 



(21) 



The joint density of states between the CB host matrix 
components and valence band vi can then be obtained 
using 



-- 1 



Gkk + Evi) 



(22) 



where Eyi is the energy of the valence band vi which, as in 
the previous section, can be the LH, HH, or split-off (SO) 
band. We take into account the nonparabolicity of the LH 
band given by Equation 13, but assume parabolic heavy- 
hole and split-off bands, as we find that the parabolic 
split-off band dispersion in the relevant energy range is 
very close to that obtained when nonparabolicity effects 
are also included. 

Having the joint DOS the optical absorption spectrum 
can be calculated using 

^>^tot = Fex ^ ^ (/c,LH + /c,HH + /c,SO ) ♦ (23) 

HrCeomQ nco 

The Greens function given by Equation 19 ignores the 
nonparabolicity of the host semiconductor CB. In order 
to consider the Kane non-parabolicity, the Green's func- 
tion given by Equation 19 is modified to include the 
conduction-valence band interaction 



Gu(£) = 



E-Ec 



E-Ey 



E 



E-Ej- AEj(E) 



-1 



(24) 



Results and discussion 

Here, we present the absorption spectra calculated using 
the five-level BAC and LCINS model and compare them 
with experiments. Perlin et al. [5,6,22-24] measured the 
optical absorption spectra for Ino.04Gao.96No.01Aso.99 and 
Ino.08Gao.92No.015Aso.985 and compared them with GaAs 
absorption data. Turcotte et al. [16,25] recently mea- 
sured the optical absorption spectrum of GaN^^Asi.^^ and 



liiyGdii.yNxAsi.x for several values of ;v and 3/. Here, we cal- 
culate the absorption spectra for Ino.04Gao.96No.01Aso.99 
and compare them with Skierbiszewski measurements at 
different temperatures. 

Five-level model for InyGai-yNxAsi-x 

The interaction between the InGaAs valence and con- 
duction bands and isolated and pair N states in 
lUyGs-i-yNxAsi^x can be described using Equation 9. The 
band structure parameters for In^Gai-^As are taken to 
vary with In composition, y, and temperature, T, as shown 
in Tables 1 and 3. Also, the energy and the interaction of 
isolated and pair N states are taken to vary with In com- 
position as given in Table 1. Figure 1 shows the calculated 
band structure of Ino.04Gao.96No.01Aso.99 where the three 
conduction sub-bands {Eu(k), E^iik) and £/(/:)) are deter- 
mined as the eigenvalues of Equation 9. Also, we consider 
the lowest eigenvalue of Equation 9 as the split-off band 
energy (£so)- The non-parabolic light-hole (£lh) is given 
by Equation 13, and the heavy-hole band (^hh) has been 
taken to be parabolic. 

The fractional T character of the conduction sub-bands 
is also required in order to calculate the joint density of 
states between the F-like conduction band components 
and the valence bands. Figure 2 shows the calculated F 
character of the CB for Ino.04Gao.96No.01Aso.99, obtained 
using Equation 15. It is observed that in the lower sub- 
band, /r has its maximum value at the CBE and decreases 
toward zero at the top of the lowest band. It increases 
again from zero to a maximum value of around 0.4 and 
goes back to zero in the middle band. Then in the upper 
band, it increases gradually from its minimum at the bot- 
tom of the upper band, approaching an approximately 
constant value around E = 2.1 eV. 



Table 3 Electrical and optical parameters In GaAs 



Parameter 


Symbol 


Values 


Lattice constant [19] 


oo 


5.65 + 3.88 X ]0-\Te - 300) (A) 


Electron effective mass [1 3] 




300 K/10K 


Conduction 


m* 


0.063/0.067(mo) 


Liglit liole 


mi 


0.076/0.082(mo) 


Heavy liole 


mh 


0.50/0.51 (mo) 


Split-off 


rriso 


0.1 45/0.1 54(mo) 


Energy gap [T = 0) 


EgO 


1.519 eV 


Varshni parameters [19] 






(Eg -Ego ^p^+Te)) 


OLj 


5.408 X 10-4 (eV/K) 






204 


SO splitting energy [19] 


Aso 


0.34 eV 


Refractive index [13] 




3.255(1 +4.5 X lO-^Te) 


Static dielectric constant [26] 


K 


12.4(1 + 1.2 X 10-^7^) 
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E(eV) 

Figure 2 The fractional T ciiaracter for InyGai-yNxAsi-x with 

y = 4% and x = 1 %. The fractional r character calculated using the 
five-level BAC model for conduction sub-bands, at 7 = 1 0 K. 



Figure 3 shows the calculated contributions of the dif- 
ferent transitions to the total absorption spectrum. The 
solid, dashed and dotted lines in this figure represent the 
contributions due to transitions from the LH, HH and 
SO bands, respectively, to the conduction sub-bands. The 
red, blue and green lines indicate transition to the upper, 
middle and lower conduction sub-bands, respectively. The 
summation of these nine transitions is shown by the black 



dash-dotted line in this figure. Multiplying this by Fex 
gives the total absorption spectrum shown by the brown 
circles in this figure. 

The usual BAC model predicts a gap in the DOS [2,27] 
of (In)GaNAs alloys. However, it is clear from Figure 3 
that the joint DOS for different transitions overlaps and 
fills this gap. Therefore, no gap is found in the absorption 
spectrum in (In)GaNAs alloys when using the BAC model. 

Figure 4 compares the calculated absorption spectrum 
using the five-level BAC model with that measured and 
calculated by Skierbiszewski [5]. The black dots in this 
figure are the experimental results for the absorption 
spectrum of Ino.04Gao.96No.01Aso.99 at T = 10 K. The 
solid black line shows the calculated absorption coeffi- 
cient using the two-level BAC model [6] with constant 
Vmc = 2.7 eV and = 1-65 eV. This line shows some dis- 
crepancies with the experimental data, especially around 
the transition to the upper sub-band of the BAC model. 
The arrows in this figure indicate the different transitions 
from the split-off, heavy- and light-hole band edges to 
the lower and upper sub-band edges, in the BAC model. 
The dashed red line in this figure displays the absorp- 
tion spectrum calculated using the five-level BAC model, 
which shows much better agreement with the experi- 
mental measurements. This model still shows some steps 
corresponding to the transitions from the HH and LH 
band to the lower, middle and upper conduction sub- 
bands (see Figure 3), whereas the experiment shows a 




E(eV) 

Figure 3 The absorption spectrum for InyGai-yNxAsi-x/ with 
y = 4% and x = 1% calculated using the five-level model at 
7=10 K.The contribution of the transitions between three valence 
bands and three conduction sub-bands are shown. The solid arrows 
designate the transitions from HH and LH bands to lower (I), middle 
(m) and upper (u) conduction sub-bands. The dashed arrows indicate 
the transitions from spin-orbit split-off band to conduction sub-bands. 



X 10^^ 




Figure 4 The absorption spectrum for InyGai-yNxAsi-x with 
y = 4%andx = 1%.The red dashed line and blue squares 
represent the absorption spectrum calculated using the five-level 
BAC model and the SCGF including the LCINS distribution of localised 
states (obtained forGaNo.oi2Aso.988)/ respectively. The black dots and 
line show the measured and calculated spectrum by Skierbiszewski 
[5] at temperature 10 K. 
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much smoother absorption spectrum and has only one 
pronounced step around E = 1.85 eV. This is due to 
the fact that in the five-band model, we have considered 
isolated and pair N states. Considering the full distribu- 
tion of N states makes the calculated absorption spectrum 
smoother, in better agreement with the experimental data. 

LCINS approach for IriyGai-yNxAsi-x 

In order to calculate the absorption spectrum using 
the LCINS model, we first calculate the Greens func- 
tion for the CB, given by Equation 24. The inset in 
Figure 5 shows histograms of the distribution of localised 
states for GaN^^Asi-^c with x = 0.84% and x = 1.2%, cal- 
culated using the LCINS approach [3,4]. This figure 
shows that the LCINS distributions for x = 0.84% and 
X = 1.2% are very similar. This implies that the LCINS 
distribution for GaN^^Asi-^c with x = 1.0% can be approx- 
imated by the one for x = 1.2%. However, the calculated 
CBE 2it x= 1.2% (indicated by £_) is about 30 meV 
lower than the CBE for x = 0.84%. In the five-band BAC 
model, we assumed that the energy gap in In^Gai-^As 
is given by Eg^caAs - 1-333/ + 0.27/. Also, the BAC 
model parameters in Table 1 suggest that including 4% In 
in Ino.04Gao.96No.01Aso.99 reduces the interaction Vn(E) 
by about 5%, which we ignore in the LCINS model. 
Because we do not have the LCINS distribution for 
Ino.04Gao.96No.01Aso.99, we approximate it here by the 
LCINS distribution calculated for GaNo.012Aso.988- 

We then solve Equations 20 and 24 self-consistently [2] . 
Figure 5 shows the CB DOS calculated by Equation 21 for 
Ino.04Gao.96No.01Aso.99 alloy. We observe that use of the 



LCINS distribution of states inhibits the gap predicted by 
the BAC model in the DOS of (In)GaNAs alloys. 

The blue squares in Figure 4 show the calculated absorp- 
tion spectrum at T = 10 K including the full LCINS distri- 
bution of N states, which is compared with the absorption 
spectrum measured by Skierbiszewski [5]. Clearly, the 
sharp steps that we saw in the five-level BAC model dis- 
appear due to the inclusion of the distribution of localised 
states. This gives a better overall agreement with the 
experimental data. The remaining discrepancies between 
the calculated and experimental data may be partly due to 
the fact that we have approximated the N distribution by 
the one that was previously calculated for GaNo.012Aso.988- 

The room temperature absorption coefficient calculated 
from the SCGF method including the full LCINS dis- 
tribution of N states is shown in Figure 6. The solid 
black line in this figure displays the absorption spectrum 
measured by Skierbiszewski [5] at T = 300 K. The red 
circles in this figure show the absorption spectrum calcu- 
lated in the LCINS model, where Fqx is calculated using 
Equation 5. The blue diamonds here indicate the result 
when we consider Fex = 1- The dashed blue and red 
lines in this figure show the optical absorption calculated 
using the five-level BAC model, assuming Fex = 1 and 
given by Equation 5, respectively. Here, we again observe 
that the results calculated using the LCINS model have 
lower values in comparison with those calculated by the 
five-level BAC model. This is because of the differences 
in the band non-parabolicity that we have considered for 
the valence bands, in the LCINS and five-level LCINS 
models. Figures 4 and 6 suggest that Fex might have 
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Figure 5 The density of state for InyGai-y NxAsi-x, witii y = 4% 
andx = 1%. Density of states calculated using the SCGF approach 
and LCINS distribution of N states at room temperature. Inset shows 
calculated distribution of N cluster state energies at low temperature, 
weighted by their interactions with the conduction band edge state 
for GaNxAs] _x with x = 0.84% and x = 1 .2%. 
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Figure 6 Room temperature absorption spectrum for 
lno.04Gao.96No.01Aso.99 calculated using LCINS and five-level 
models. The absorption spectrum calculated using the LCINS (red 
circles and blue diamonds) and five-level BAC (dashed lines) models. 
The red and blue curves display the results with and without 
including the electron-hole interaction. The solid black line shows the 
experimental data. 
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a stronger temperature dependence than what we have 
considered in our calculation. The temperature depen- 
dence is considered only in the static dielectric constant as 
shown in Table 3 in calculating exciton Rydberg energy in 
Equation 6. 

The absorption spectrum for GaNxSbi_x 

We can apply our method to calculate the band structure 
and absorption spectrum of other dilute nitride alloys. 
Here, we extend our calculations to investigate the absorp- 
tion spectrum of GaNSb. The room temperature band 
gap of GaSb is about 725 meV, around half that of GaAs. 
Lindsay et al. [17,28] have reported that N-related defect 
levels lie close to the CBE in GaNSb and therefore strongly 
perturb the lowest conduction states in this alloy. The 
band gap and optical properties in GaN^^Sbi-^c have been 
shown to be strongly affected and highly sensitive to the 
distribution of the nitrogen atoms. Lindsay et al. [28] 
found that there is a wide distribution of N levels lying 
close to and below the CBE. The higher- lying N states 
push the CBE down in energy, as in GaAs, but the large 
number of lower energy N states are calculated to mix in 
strongly with the conduction band edge states, severely 
disrupting the band edge dispersion in GaNSb. 

Here, we first investigate the band structure and opti- 
cal absorption spectra of GaN^^Sbi-^^ in the five-level BAC 
model and compare the results with the absorption spec- 
tra measured by Veal et al. [20] and Jefferson et al. [29]. 
We then apply the SCGF method to GaN^^Sbi-^^; in the 
Section 'LCINS model for GaN;cSbi_;c- As the LCINS dis- 
tributions have not yet been calculated for these alloys, we 
modify those calculated for GaNAs alloys and use them in 
our calculations. 

Five-level BAC model for GaNxSbi_x 

When a single Sb atom is replaced by N in GaSb, the N 
atom introduces a localised state with energy Ejsf, How- 
ever, a GaNSb alloy can also contain clusters of N atoms, 
such as N-N nearest neighbour pairs as well as larger 
clusters that introduce states in the band gap of GaSb. 
Table 2 contains the band parameters that we use for 
GaN;cSbi_;c> including the isolated N state energies, Ejsf 
and N pair state energies, £aqv relative to the valence 
band maximum energy and the BAC interaction parame- 
ters and PmN' As shown in this table, isolated N states 
are calculated to be less than 0.1 eV above the conduc- 
tion band minimum, while the N pair states have energies 
that lie in the GaSb band gap. The calculated energy 
gap of GaN;cSbi_;c depends strongly on the assumed N 
distribution, reflecting that N cluster states introduce a 
series of defect levels close to the CBE in this alloy. In 
addition, the interaction parameters {pM and Pnn) in 
GaN;cSbi_;c are calculated to be about 20% larger than for 
GaNAs alloys. 



Figure 7 shows the conduction and valence band dis- 
persion, calculated using the five-level BAC model given 
by the Hamiltonian of Equation 9. The solid lines in this 
figure show the conduction sub-bands. Here, we include 
the isolated and pair N states and their interaction with 
the GaSb conduction and valence bands, as explained in 
the previous section. Since Ejsf is very close to the GaSb 
CBE and £aqv just below it, we observe that the lower sub- 
band (£/) is almost flat and located within the GaSb band 
gap. The band edge minimum for this band is 0.39 eV and 
its maximum energy is 0.45 eV. This implies that substitu- 
tion of only X = 1.2% N by Sb in GaSb reduced the energy 
gap rapidly from 725 to 390 meV. This value for the band 
gap of GaN;cSbi_;c with x = 1.2% is very close to that 
which was previously measured [20] and calculated using 
k.p [28] and ab initio pseudopotential [30] calculations. 

The middle sub-band lies between 0.55 and 0.78 eV, 
and the upper sub-band {Eu) minimum is close to 1.0 eV. 
The blue dashed line in Figure 7 shows the non-parabolic 
spin-orbit split-off band (£so)> calculated by the lowest 
eigenvalue of Equation 9. Also, the non-parabolicity of the 
light-hole band (£lh) has been taken into account using 
Equation 13, while we assumed that the heavy-hole (£hh) 
band has a parabolic dispersion. 
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Figure 7 The band dispersion of GaNo.012Sbo.988 calculated by 
the five-level BAC model. The solid lines display the conduction 
bands including upper {Eu), middle {Em) and lower {Ej) sub-bands. The 
dashed lines show the spin-orbit split-off {Esq), light-hole (^lh) and 
heavy-hole (^hh) bands. 
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Given the band dispersion, we can calculate the opti- 
cal absorption as described earlier for InGaNAs. The red 
dashed line (with circles) in Figure 8 shows the absorption 
spectrum of GaNo.012Sbo.988 calculated using the five-level 
BAG model. Green and brown solid lines in this figure 
show the absorption spectra measured by Mudd et al. 
[31] for GaN;cSbi_^ with x = 1.18% and x = 1.22%. Our 
calculated absorption edge is in good agreement with 
these experiments. However, there are two sharp steps in 
the calculated spectra corresponding to transitions from 
the light- and heavy-hole bands to the middle and upper 
sub-bands. We observe that the experimental absorption 
spectrum, a(E), increases from zero at the CBE, to about 
3 X 10^ cm~^ at energy E = 0.55 eV, that is the calcu- 
lated band edge of the middle sub-band. After this point 
the slope of the spectra decreases up to £ = 1 eV, the 
minimum of the upper conduction sub-band. Then, due 
to the transition from valence bands to the upper conduc- 
tion sub-band, the magnitude of the calculated absorption 
spectra increases rapidly. The spin-orbit splitting energy, 
Aso, is 0.76 eV. Therefore, transitions from the split-off 
band commence at 1.15 eV, where we see a small increase 
in the calculated absorption spectrum due to transitions 
from the spin-orbit split off band to the lowest conduc- 
tion sub-band E\, The calculated optical absorption using 
a 5 X 5 k.p Hamiltonian, accounts well for the absorption 
edge. Wang et al. [32] have also measured the absorption 
edge of GaN;cSbi_;^; with x = 0.3%, x = 0.7% and x = 
1.4%, with the measured band edge energies in very good 
agreement with those calculated by the five-level model of 
this paper. 
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Figure 8 The absorption spectrum of GaNo.012Sbo.988- The 

absorption spectrum for GanxSb] _x with x = 1 .2%, calculated using 
the SCGF model and including the distribution of localised states 
(blue diamonds), in comparison with experimental data measured for 
GaNxSbi_x withx= 1.18%andx= 1 .22% (solid lines). The results 
calculated by the five-level BAG model are displayed by the red line 
with circles. 



We note however that there are two sharp features in the 
calculated results that are not observed in experiments. 
This could be due to the fact that we included only isolated 
and pair N states in this model and ignored the distri- 
bution of N states and their inhomogeneous broadening. 
The results for the calculated absorption spectrum using 
a five-level BAC model suggest that we need to include 
the full distribution of N states in optical absorption 
calculations. 

LCINS model for GaN^Sbi _x 

It has been shown that the calculated electronic structure 
of GaN;cSbi_;c strongly depends on the assumed distribu- 
tion of N atoms [28] . Therefore, in order to calculate an 
accurate band dispersion for this alloy, we need to have 
the distribution of localised states. Unfortunately, such 
a distribution has not been calculated for GaN^^Sbi-^c- 
However, we expect that the distribution of N states in 
GaN;cSbi_;c should have a similar general form to the 
LCINS distribution that Lindsay et al. [3] have calculated 
for GaN;cAsi_;c and for GaN;^^-;^ alloys [17]. Therefore, 
here, we consider the LCINS distribution of N states in 
GaN;cAsi_;c and, with some small modifications, use that 
for GaN;cSbi_;c alloys. 

In the previous section, we have seen that the cal- 
culated energy of an isolated N state is at about 
0.82 eV. So, we first need to shift the LCINS distribution of 
GaNo.012Aso.988 to locate the highest peak at this energy. 
The dashed red line in the inset of Figure 9 displays the 
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Figure 9 The DOS of GaNo.012Sbo.988/ calculated using the SCGF 
method and including the LCINS distribution given in the inset. 

Inset displays the distribution of N states assumed for GaNo.012Sbo.988 
(blue line), in comparison with the LCINS distribution of GaNo.012 
Aso.988/ shifted down in energy by 888 meV (dashed red line). The 
zero of energy is taken to be at the top of the GaSb valence band. 
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LCINS distribution of N states, weighted by V?/^s/^> cal- 
culated for GaN^Asi_;c with x = 1.2%, and shifted down 
in energy by 888 meV. This distribution can be approxi- 
mated by three Gaussian distributions, each correspond- 
ing to different N environments. It is observed in this 
figure that if we align the main peak at £ = 0.82 eV, the 
lowest peak corresponding to pair N-N states is located at 
0.55 eV, which is higher than the values that we consid- 
ered for Emn in the BAG model. Therefore, we shift the 
Gaussian distributions corresponding to pairs and larger 
clusters of N states down by a further 70 meV. More- 
over, the BAG model parameters in Table 2 suggest that 
in GaNSb, the interaction parameters, and Pnn> are 
20% stronger than in GaNAs. Therefore, we multiply the 
N LGINS values by 1.44 to account for this difference. The 
blue solid line in the inset of Figure 9 presents the distri- 
bution of N states that we consider for GaNo.012Sbo.988 in 
our calculation. 

Having the distribution of N states, we are able to 
calculate the Greens function for GaNo.012Sbo.988 using 
Equations 19 and 20, self-consistently. Also, the density of 
CB states can be calculated using Equation 21. Figure 9 
shows the DOS of GaNo.012Sbo.988 calculated by the SCGF 
method and including the distribution of localised states 
shown by the solid blue line in the inset of Figure 9. The 
gaps corresponding to isolated and pair N states are clearly 
observed in this plot. Also, at energies around 0.65 eV, the 
DOS has a small gap that is related to the higher cluster of 
N states. 

We can also calculate the absorption spectrum using the 
SCGF model. The blue line with diamonds in Figure 8 
shows the calculated absorption coefficient using this 
method. As expected, this method shows a better agree- 
ment with experiments than the result of the five-level 
BAG model (shown by red circles in this plot). 

For the considered N distribution, this calculation sug- 
gests more gaps in the DOS of GaNo.012Sbo.988 compared 
to GaNo.012Aso.988 [2]. However, experimental data indi- 
cate that there are fewer features in the GaNSb absorption 
spectra than in the GaNAs ones. This could be due to 
inhomogeneities in the samples investigated experimen- 
tally, either due to fluctuations in the N composition in 
the experimental samples or because of intrinsic differ- 
ences between the short-range N ordering in GaNSb and 
in InGaNAs samples. 

Recent work by Mudd et al. [31] has shown that the 
composition dependence of the energy gap in GaN^^Sbi-^^ 
is well described using a three-level model including inter- 
actions between the host matrix band edge and the N 
isolated states and N-N pair states. The energy gap cal- 
culated using the LCINS model is also determined pri- 
marily by these interactions. The energy gap calculated 
here using the SCGF and LCINS method is consistent 
with experiment for the N composition x = 1.2% which 



we consider and should closely follow the theoretical 
energy gap results presented in [31] as a function of N 
composition x. 

Conclusions 

In this paper, we presented an analysis of the optical 
absorption spectra of dilute nitride alloys, calculated using 
the band structure model presented in our earlier work 
[2]. We have considered two different models to calcu- 
late the absorption spectra in InGaNAs and GaNSb alloys 
and compared our results with experimental measure- 
ments. We note however that there are some discrepan- 
cies between experimental data in similar samples that 
make quantitative comparison difficult. 

Two models have been considered to calculate the 
absorption spectrum in these materials: a five-level BAG 
model and a LCINS -based model. The five-level BAG 
model included isolated and pair N states and their inter- 
actions with the host semiconductor valence and con- 
duction bands. The results of this model for InGaNAs 
alloys give an overall good agreement with experiments, 
and predict accurate absorption edge for these alloys. 
However, the results of the five-level BAC model include 
several additional features not seen experimentally, sup- 
porting the need to consider a full distribution of N state 
energies in the electronic structure calculations. 

We therefore extended our calculations to include 
the LCINS distribution using the SCGF approach pre- 
sented in [2]. The calculated absorption spectra using 
this approach for InGaNAs provide very good agree- 
ment with experiments, supporting the validity of the 
LCINS approach to describe dilute nitride conduction 
band structure. 

Our calculated absorption spectra for GaNSb alloys fit 
well with experiments at the absorption edge [31], and 
predict the correct band gap in these alloys. However, 
the absorption spectrum calculated in the BAC model 
contains features associated with individual transitions to 
lower and upper sub-bands in the model that are not seen 
in the measured absorption spectra. Taking the distribu- 
tion of localised states into account reduces the impact 
of these features and gives results more similar to exper- 
imental absorption. But we still see some dips in our 
calculated spectra that are not seen in any experiment. We 
conclude that the distribution of N states in the GaNSb 
alloys studied are different from that for InGaNAs sam- 
ples. We conclude that further work is required to address 
and resolve why more structure is found in the calculated 
absorption spectra compared to what is observed in the 
experimentally measured spectra. 
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